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Abstract 



In this paper we investigate the problem of searching monotone multi-dimensional 
arrays. We generalize Linial and Saks' search algorithm [2] for monotone 3-dimensional 
arrays to d-dimensions with d > 4. Our new search algorithm is asymptotically op- 
timal for d = 4. 
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1 Introduction 

In this paper, we investigate the problem of searching monotone multi-dimensional 
arrays. Suppose we are given a d- dimensional array with n entries along each 
dimension 

A n ,d = { a ii,i2,-,*dl*i> *2, • • • , id — 1) 2, . . . , nj. 

We say that the array A n ^ is monotone if its entries are real numbers that 
are increasing along each dimension. More precisely, if i\ < j±, i 2 < ji,- ■ ■ , 
id < 3d then Oij,^,...,^ < a h,n,...,j d - 111 other words, if P = [n] d is the product 
of d totally ordered sets {1,2, ... ,n}, than A n ^ d is consistent with the partial 
order P. 

The search problem is to decide whether a given real number x belongs to 
the array A n ^ by comparing x with a subset of the entries in the array. The 
complexity of this problem, denoted by r(n, d), is denned to be the minimum 
over all search algorithms for A n ^ of the number of comparisons needed in the 
worst case. Note that for d — 1, this problem reduces to searching a totally 
ordered set. In this case, the binary search algorithm is optimal and requires 
at most |Tog 2 (n + 1)] comparisons in the worst case. 

We first briefly review some previous work. In [3], Linial and Saks presented 
some general results on the complexity of the above class of search problems, 
for any finite partially ordered set P. In [2] they studied the problems for 
general finite partially ordered set P and also gave more precise results for the 
case where P = [n] d , for dimensions d > 2. They observed that for d — 2, it 
had been known that r(n, 2) = 2n — 1 [1]. For the generalized case d > 2, they 
showed that the order of r(n, d) is 0{n d ~ 1 ). More specifically, they proved that 
for d > 2, 

ci(d)n d_1 > r{n,d) > c 2 (d)n d ~ 1 + o{n d " 1 ) , 

where c\{d) is a nonincreasing function of d and upper bounded by 2, and 
c 2 (d) = J(24/Tr)d" 1 / 2 + o(d- 1 / 2 ). The upper bound ci(d)n d_1 was obtained by 
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using a straightforward search algorithm which partitions A n ^ into n isomor- 
phic copies of A Ut d-i and searches each copy separately. They also described a 
more efficient algorithm for d = 3 and proved the following bounds on r(n, 3): 

,3n 2 , . . 3n 2 

< T \ n i 3) < — — h cn In n. 

In the above inequality, c is a positive constant, and so the bounds are asymp- 
totically tight. An open problem left is whether their search algorithm for 
d = 3 can be generalized to higher dimensions. 

In this paper, we present new search algorithms for monotone <i-dimensional 
arrays for d > 4, by generalizing the techniques in [2] to higher dimensions. 
For d > 4, the search complexity of our algorithms is 

r(n,d) < -^—n d - 1 + 0(n d - 2 ). 
d — 1 

The above bound is tight for d — 4, up to the lower order terms. 

The rest of the paper is devoted to the description and analysis of the new 
algorithms. We start with the case where d — 4. This special case best il- 
lustrates the main idea, and it is also easier to visualize the subspaces that 
are encountered in the search algorithm. Then we describe the generalized 
algorithm for d > 4. 

Before presenting the technical details, we describe some basic notation and 
convention that we will follow throughout the paper. In general, we use capital 
letters to represent sets and small letters to represent numbers. The sets that 
we need to consider are often subsets of A n ^ for which some of the subscripts 
are fixed, and we use some simple notation to represent them. For example, 
we use Q = {o>i,i2,i 3 ,u} to denote a "surface" of the 4-dimensional array A n 4 
for which the first subscript of a is fixed to be 1. It is understood that all other 
subscripts range between [1, n], and we often omit the specification "i 2 , i 3 , i 4 = 
1, 2, . . . , n" if it is clear from the context. 



2 Searching 4-Dimensional Arrays 

In this section, we present a |n 3 + 0(n 2 ) algorithm for searching monotone 
4-dimensional arrays. The algorithm is optimal up to the lower order terms. 

We start with a lower bound on r(n, 4) which will be seen asymptotically tight 
later, followed by the description of an algorithm for partitioning monotone 
two-dimensional arrays, which will be a useful subroutine for our searching 
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algorithm. Then, we will present the main idea and the details of our search 
algorithm for 4-dimensional arrays. 

2.1 A lower bound on r(n, 4) 

Using the method in [2], we can calculate a lower bound on r(n, 4). Let [n] 
denote the totally ordered set {1,2,..., n}, and let 

Z>i(ra, 4) = {(ij, i 2 , i 3 , u) G [n] 4 |«i + i 2 + i 3 + i 4 = 2n + 1}, 
D 2 (n, 4) = {(ii, i 2 , i 3 , u) e [n] 4 |«i + i 2 + «3 + «4 = 2n + 2}. 

Define D(n,4) = D 1 {n,4)UD 2 {n,4). Then £>(ra,4) is a section (see [2]) of [n] 4 , 
and there is no ordered chain having length more than 2 in D(n, 4). Therefore, 
r(n, 4) is lower bounded by \D{n, 4)|. Let 

= {(ii, «2, «3, «4) £ [2ra + 1] 4 | «i + ii + *3 + ^4 = 2n + 1}, 
= {(ii, «2, «3, «4) G X| i fc > n} for fc = 1, 2, 3, 4, 

Z = {(h, *2, «3, G [™ + 1] 4 | *1 + *2 + ^3 + U = n + 1}. 

It is easy to see that |Yjt| = |Z| = (^j for fc = 1,2,3,4. Thus, \Di(n,4)\ = 
\X\-ELi \Y k \ = (?) -4g) = §(2n 3 -2n). Similarly, |L> 2 (n,4)| = ( 2 " 3 +1 ) - 
4( n + x ) = |(2n 3 + n). Therefore, 
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r(n,4) > \D(n,4)\ = |£>i(M)| + |£> 2 (n,4)| = -n 3 - -. 

Partitioning 2- dimensional arrays 

In [2], Linial and Saks gave a simple search algorithm for an m x n matrix 
(m,n > 1) with entries increasing along each row and column. The algorithm 
needs at most m + n — 1 comparisons. We will refer to this algorithm as the 
Matrix Search Algorithm. Since A nj2 is isomorphic to an n x n matrix, we can 
adapt the Matrix Search Algorithm to partition A n 2 into two subsets S and L 
given an input x, such that S contains entries smaller than x and L contains 
entries larger than x, using at most 2n — 1 comparisons. Below, we provide the 
detailed description of the partition algorithm for the sake of completeness. 

Algorithm: Partition 2-Dimensional Array 

Input. 
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Fig. 1. u,v,S,L: partition of the monotone 2-dimensional array A n ^ 2 

• A real number x. 

• A monotone 2-dimensional array A n>2 = {ou,j 2 }- 
Output. 

• If x G Ai,2, output (h,^) such that = x. 

• If x £ A n>2 , output a partition {u, v, S, L} of A Uj2 with the following prop- 
erties: 

- u and v are two arrays each contains n integers such that i\ < u[i 2 ] iff 
Oil ,i 2 < ^ an d ^2 < iff aii,j 2 < x - 

- S and L form a partition of «2)|*i, «2 £ N} such that if (ii, 22) £ 5* 
then aj l5 j 2 < x, and if (ii , 22) £ £ then > x. 

Procedure. 

• Initially set S = L = 0. 

• View A n>2 as an n x n matrix and repeat comparing a; with the element e 
at the top right corner of the current matrix. 

- If x > e, then eliminate the first row of the current matrix and put their 
entries into S; 

- If x < e, then eliminate the last column of the current matrix and put 
their entries into L; 

- If x — e, then return this entry and exit. 

• Stop when the partition is finished, thus also obtain u and v (see Fig. 1). 

We will use the notation u,v,S,L throughout the paper. Sometimes we will 
introduce subscripts to them to represent the dimension indices to be con- 
sidered. Ignoring the indices, these four variables have the following useful 
relations: 



S = {(ii,i 2 )|l < k < u[i 2 ]} = {(i 1 ,i 2 )\l < i 2 < v[h]} 



L = {(ii,i 2 )\u[i 2 ] < h < n} = i 2 )\v[h] < i 2 < n}. 



(2) 
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Obviously, S fl L = <fi and S U L = [n] 2 , hence \S\ + \L\ = n 2 . In addition, 
\S\ = u[l] + . .. + u[n] = v[l] + . ..+v[n] and \L\ = (n-u[l]) + . .. + (n-u[n\) = 
(n — v[l\) + . . . + (n — v[n\). 

Notice that when m = or n = 0, we can "search" an m x n matrix using 
comparisons. Therefore, based on the Matrix Search Algorithm in [2], we have 
the following lemma that will be useful later. 

Lemma 2.1 For m,n > 0, any m x n matrix with entries increasing along 
each dimension can be searched using at most m + n comparisons. 

Proof. If m = or n = 0, the matrix is empty, thus needs no comparison. If 
m, n > 0, using the Matrix Search Algorithm, we can search the matrix using 
at most m + n — 1 comparisons. Therefore, the lemma holds. ■ 



2.3 Main idea of the search algorithm 

The main idea of our algorithm for d = 4 is to first search the "surfaces" 
(3-dimensional arrays) of A n 4 and then the problem reduces to searching a 
"smaller" array A n _ 2 ,4- At a high level, searching the surfaces consists of two 
major steps: 

• Step 1: Select 8 special 2-dimensional arrays and partition each into two 
subsets L and S, where elements in L are larger than or equal to x, and ele- 
ments in S are smaller than x, using the algorithm Partition 2-Dimensional 
Array. 

• Step 2: Search the 8 "surfaces" of A n> ±. The subsets S, L obtained in Step 
1 help to "cut" each surface into a sequence of 2-dimensional matrices that 
allows searching with less comparisons. 

2.4 Description and analysis of the search algorithm 

Now we are ready to present our search algorithm for d — 4. As explained 
in Section 2.3, the algorithm is recursive, which reduces n by two for each 
recursion. Without loss of generality, we consider the case where x ^ A n ^. 
We first describe the algorithm and then analyze the number of comparisons 
needed. 

Step 1. Apply the algorithm Partitioning 2-Dimensional Array to divide each 
of the following eight 2-dimensional arrays into two subsets (the eight arrays 
are defined by fixing two of the subscripts to either 1 or n, thus reducing the 
number of dimensions by two). 



6 



M l = {o»i,ia,l,n} 
M2 = {ttn,i 2 ,i3,l} 

M 4 = {0^,1^,14} 



5*1, L\\ 
S2, L2; 



M l = { a h,i2,n,l} 

M% — { a l,i2,i:i,n} 

= {a n ,l,i 3 ,i 4 } 

m: 



Q* J * 
Q* T * . 

<-> 2 , Ll 2i 

S3, LI; 



Q* T * 
4 , lv 4 . 



The eight pairs of "mutually complementary" subsets S^, L fc and S^*, L k * (k 
1, 2, 3, 4) have the following properties: 



a ii,?2,l,n < < a jiJ2,l,r. 



for (ii,i 2 ) G 5i and (ji,j 2 ) G £1 



,j2j3,i f° r (^2^3) G S 2 and {j2ijz) G -^2 

ai,n,i 3 ,i4 < x < a i,n,i3,i4 for (*3,«4) € S3 and (j 3 , j 4 ) G L 3 



^ii,l,n,i4 < < ^?i,l,n 



for (i 4 ,ii) G S 4 and (j 4 , Ji) G L A 



a 



,i a ,n,i < » < ajua.n.i for (ii,^) G and (ji, j 2 ) G L* 



ai,i 2 ,i 3 ,n <x < aij 2 j3,n for (12,13) G and (j 2 ,j 3 ) G L* 

a n, 1,13,14 < 3; < a "ri,l,j: i ,j4 

for (i 3 ,i 4 ) G S3* and (j 3 , j 4 ) G L* 

a h,n,l,u < x < a ji,n,l,j4 

for (i 4 ,ii) G S 4 and (j 4 , Ji) G £ 4 

In addition to the eight pairs of subsets, the algorithm also outputs u^, and 
ul,v%., corresponding to Sk,Lk and Sk*,Lk* respectively, with the properties 
given in Equation 1 and 2. For each k, at most 2n — 1 comparisons are needed 
to partition M k (M£). Thus, at most 8 x (2n — 1) comparisons are needed in 
this step. 

Step 2. Search the following eight 3-dimensional surfaces of A n ^ (each surface 
is defined by setting one of the subscripts to either 1 or n, thus reducing the 
number of dimensions by one). 

Ql = { a l,«2,«3,«4} Ql = { a n,i 2 ,«3,«4} 

Q2 — {°il, 1^3,14} Q2 = { a «l,™,«3,«4} 

Q3 = {^11,12,1,14} Q3 = { a u,«2,n,«4} 

Qi — { a il,«2,«3,l} Qi = {°n,«2,«3,™} 

By symmetry, we only need to show how to search Q\. The algorithm proceeds 
by fixing i 3 = i 3 for i 3 = l,2,...,n and searching each of the 2-dimensional 
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(a) 



(b) 



(c) 



Fig. 2. Searching the 3-elimensional surface Q\ = {ai,i 2 ,i 3 ,i 4 } of A n ^. (a): partition 
of M| = {ai 5 j 2) i 3)n } into S% and L 2 ; (b): partition of M3 = {ai ;n< i 3: i 4 } into 53 and L3; 
(c): the 'pyramid' composed of a sequence of 2-dimensional matrices to be searched. 

array {a li2i ' i4 }- A useful observation is that for each i' 3 , we can restrict the 
search to a smaller matrix (in contrast to an n x n matrix) by leveraging on 
information obtained in Step 1. 



Below, we explain the above observation and Step 2 in more details. Consider 
an element a, ■ ■> ■ G Q\. If (i2,H) G SZ, then we know that a, , / . < 

1,12,1^,14, ^ 1 \ z ' 0/ z ' 1,12,«3,«4 — 

Oi / „ < Hence, in order for a, . / • = x, it must be the case that (12, u) G 

L 2 , or equivalently, ^2^3] < 12 < n - Similarly, we can conclude that in order 
for a l i2 j i4 = x, it must be the case that (i 3 , i 4 ) G L 3 , or equivalently, 1*3 [i 3 ] < 
24 < n. Hence, we obtain a constraint on the indices (12, h)- By lemma 2.1, 
searching this restricted (n — u 2 [i 3 ] ) x (n — f 3 [i 3 ] ) matrix needs at most (n — 
u 2 lh] ) + ( n ~~ ^3 ih] ) comparisons. Notice that n — w^fe] is the number of entries 
(«2,*3)'s in with i 3 = i 3 , and n — ^3^3] is the number of entries (13, i4)'s in 
L 3 with i 3 = i 3 . Thus, 

(n-«2[4]) + (^-^[4]) = |{(*2,»s) e ^2^3 = 4)1 + |{(*3,i4) e £3^3 = h}- 

When i 3 ranges over 1,2, ... ,n, we obtain that the total number of compar- 
isons needed to search Q x is at most N(Qi) = \L 2 \ + \L 3 \ (see Fig. 2). 

Similarly, if a nj j 2j j 3) j 4 G Q\ equals to x, it must be the case that ^2,23) G S 2 
and (*3, 24) G S 3 , it follows that the total number of comparisons needed to 
search Q\ is at most N{Q\) = \S 2 \ + \S 3 \. 

Using similar arguments, the numbers of comparisons needed for searching the 
above eight 3- dimensional arrays are: 



N(Qi) 


= \Ls\ 


+ \L* 2 \, 


N(Ql) 


= \S2\ 


+ 1^3*1 


N(Q 2 ) 


= 




N(Q* 2 ) 


= \s,\ 


+ 1^4*1 


N(Qs) 


= \Li\ 


+ \LH 


N(Q* 3 ) 


= LS4 


+ 


N(Q 4 ) 


= 


+ \Ll\, 


N{Ql) 


= \sl\ 


+ 1^1 
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Therefore, the total number of comparisons needed for searching these eight 
arrays is at most 

4 

J2(\S k \ + \L k \ + \S* k \ + \L%\) = 4 x 2n 2 = 8n 2 . 

k=l 



Steps 1 and 2 leave an (n — 2) 4 array 

Ai-2,4 = { a ii,i2,i3,i4rl) *2, ^3^4 = 2, . . . , 71 — 1}. 

Hence, we have for n > 2, 

r(n, 4) < r(n - 2, 4) + 8n 2 + 8(2n - 1) 
From this recursion we can get (see Equation 4 for the derivation) 

r(n,4) < ii 3 + 0(n 2 ). (3) 



3 Searching d-Dimensional Arrays 

The algorithm for 4-dimensional arrays can be generalized to higher dimen- 
sions (d > 4). The main idea is quite similar: the 2d "surfaces" ((d — 1)- 
dimensional arrays) of A n ^ can be searched using 2dn d ~ 2 + 0(n d ~ 3 ) compar- 
isons. We achieve this in two steps. First, select 2d special (d — 2)-dimensional 
arrays and partition each of them into two subsets S and L. Second, we search 
the 2d "surfaces". The subsets {S, L} will help cut some part of each surface, 
i.e., reduce the comparison number. In particular, if we fix (d — 3) subscripts, 
the resulting part is a smaller matrix (in contrast to an n x n matrix). An 
a x b matrix can be searched using at most a + b comparisons (Lemma 2.1), 
adding them up for all the 2d "surfaces" , we can get the desired upper bound. 

First we describe how to select and partition the (d — 2)-dimensional arrays. 
Define Mi = {0,^,12,.. .,i d £ A nt d\id-i — 1, id = n}. Consider the case where 
x ^ Mi. For fixed i 2 = i 2 ,i3 = i 3 ,...,id-3 = ^-3 (where i 2 , i 3 , . . . , i d _ 3 E 
{1,2, . . . ,n} are constants) we can get two integer arrays u[n] and v[n] such 
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that 



*-l«=4 <«i-3=ii_3; ia-i=M«=« < X fOT fl " M[ ^ 2]; 



«n. 2 ,....J J2= 4,...,, d _3 = ^ 3 ; id _ 1= i,i d =„ > x for ix > ufid-a]; 

«H^,..^l i2 =4,..., id _3= i '_ ,; i d -l=Md=n < * fOT *d"2 < 



°ii,i2,— )*dU 2 =4,— ,id-3=i'j-v 



l=Md=n >x fori d _ 2 >u[ii]; 



by using the algorithm Partitioning 2-Dimensional Array, at most 2n — 1 
comparisons are needed for each fixed i 2 , . . . ids- Thus using at most n d ~ 4 (2n— 
1) comparisons we can get two integer arrays u\ and v\ of sizes n d ~ 3 such that 

- If ix < Mi[? 2 , • • • ,^-2], then a ilii2v .. !id | id _ 1=:Md=ri < x. 
Otherwise, a ilM ^ M \ id _ 1=1M=n > x. 

- If id-2 < • • • ,^-3], then ai 1 ,i 2 ,...,i d |i d _ 1 =i 1 i d = n < 
Otherwise, Oi^ij, i d ]i d _ 1 =i,< d =n > x. 

Thus, we can partition [n] d ~ 2 into two subsets Si and L\ such that 

~ a h,i2,...,i d \i d - 1 =l,i d =n < ^ f° r • • • 5 *d-2) G 5*1 . 
~ a ii,i2,...,idlid-i=l,»d=n > ^ ^ OT (*1> • • • ' G -^1- 

Obviously, we have 

- (ix, . . . ,i d _ 2 ) G Si if and only if ^ < m[i 2 , ■ ■ ■ ,^-2] (also i d _ 2 < . . . ,^- 3 ])- 

- (ix, . . . ,i d _ 2 ) E L x if and only if ?! > ui[i 2 , ■ ■ ■ ,^-2] (alsoi d _ 2 > • • • ,id-s])- 

Next we describe the algorithm for searching <i-dimensional arrays A n ^, for 
d > 4. Without loss of generality, we consider the case where x G" 

Step i. Partition each of the following 2d (d — 2)-dimensional arrays into two 
subsets. 

M k = {a>i 1 ,i 2 ,...,i d \ik-2 = Mfc-i — n } : Sk,L k 
M k = {aii,i2,..,i d \ik-2 = n,i k -i = 1} : S* k ,L* k 



k = 1, 2, . . . , d (here i k _ 2 means i (fc _ 2 ) mo d d, and i fe _i means i (fc _i) mod d ). 

We get 2(i pairs of mutually complementary subsets S k , L k and S k , LI with 
the following properties: 
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«ii,...,i d |i fc _2=l,ifc_i=n x a ji,-..,j d I jk-2=l ,jk-i=n 

for . . . , G S'fe and (j k , . . .Jk+d-s) G L fe ; 

a il,...,i d |i fc _2=n,i fc _i=l < ^ < a ii,...,idlifc_2=nj'fe-i=l 

for (i fc , . . . , i fe+d _ 3 ) e SI and . . . , jk+d-3) e 
fc = 1, 2, . . . , d (here i fc+(i - 3 means i(k + d_ 3 ) mod d etc.). 

For the pair S k and we have two (d — 3)-dimensional arrays u k and v k such 
that, if i k < u k [ik+i, • • • , ik+d-z] then . . . , i k +d-z) G Sfc else (i fc , . . . , ? fc+d _ 3 ) e 
L fc ; if i fc+ d-3 < ^fc[4, • • • , then (i k , . . . , i k+d - 3 ) G S k else (i fe , . . . , i fc+d _ 3 ) G 

Lfc, = 1, 2, . . . , d. Similarly, we have u* k and vl for the pair S% and L£, for 
fc = 1,2, 

In this step, we obtain 4d (d — 3) -dimensional arrays u k , v k , u* k , v k (k = 
1, 2, . . . , d), using at most 2dn d ~ 4 (2n — 1) comparisons. 

Step 2. Search the following 2d (d — l)-dimensional surfaces of A n>d , which 
are defined by fixing one of the subscripts to either 1 or n. 

Qk = {a iu ..., id \i k = 1} Q* k = {ai u ...,i d \ik = n} 

k = 1,2, 

By symmetry, we only need to consider searching Q\ = {ai,i 2 ,...,i d }- 

If a i,...,i d ^ Qi equals to x, we have (i 3 , . . . ,id) G L 3 and (i 2 , ■ ■ ■ ,i d -i) G L\. 
For fixed i 3 = i 3 , . . . , = (where i 3 ,. . . ,i d _ 1 G {1,2,..., n}), there exist 
two integers u = u^i^ ■ ■ ■ , and v = v 3 [i 3 , . . . , i d _i\ such that only when 
u<i 2 <n, (i 2 , i' 3 , . . . , i d _i) G L2; only when v < i d < n, (i' 3 , i d ^i d ) G L 3 . 
Thus for fixed i 3 = i 3 , . . . , i^-i = i^-i, on ly when u < i 2 < n and v < i d < n, 
a i i 2 i' i' i d P oss ibly equal to x. Searching this (n — u) x (n — v) matrix 
needs at most (n — u) + (n — v) comparisons. Notice that n — u is the number 
of elements (i 2 , . . . , id-i)'s in L\ with i 3 = i 3 ,...,i d _i = i d _ 1 , and n — v 
is the number of elements (i 3 , . . . , i d )'s in L 3 with i 3 = « 3 , . . . , i d _i = i d -i- 
Thus (n - u) + (n - v) = \{(i 2 , . . .,i d 1) G L*|i 3 = i' 3 , . . . ,i d _! = z^_ x }| + 
|{(i 3 , . . .,i d ) G L 3 |i 3 = i' 3 , . . . ,? d _i = irf !> I - When (i 3 , . . . ,i d _{) ranges over 
all elements in [n] a!_3 , we obtain that the total number of comparisons needed 
to search Qi is at most N(Q 1 ) = |L 3 | + \L^\. 

Similarly for Q\ = {a nii2i ...^}, we have N(Ql) = \S 2 \ + \S%\. 

Using similar arguments, the numbers of comparisons needed for searching the 
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above 2d surfaces are 

N(Q k ) = \L k+2 \ + |L* +1 | N(Ql) = \S k+1 \ + \S* k+2 \ 

k = 1, 2, . . . , d (here L k+2 means L {k+2) mod d etc.). 
Thus, searching these 2d (d — l)-dimensional surfaces needs at most 

j^{\L k+2 \ + \L* k+1 \ + 15*+! | + \S* k+2 \} = dx 2n d ~ 2 = 2dn d ' 2 

k=l 

comparisons. 

Steps 1 and 2 leave an (n — 2) d d-dimensional array 

A n - 2 ,d = {a>h,...,i d \ii, ■ ■ ■ , U = 2, . . . , n - 1}. 
Hence the generalized recursion is 

r(l,d) = l 
r(2,d)<2 d 

r(n, d) < r(n - 2, d) + 2dn d ~ 2 + 2dn d " 4 (2n - 1) for n > 2 

From the recursion, for fixed d there exist a constant C and e G {1,2} such 
that 

r(n, d) < r(n - 2, d) + 2dn d " 2 + 4dn d " 3 

< r(n - 4, d) + 2d(n d " 2 + (n - 2) d ~ 2 ) + 4d(n d ~ 3 + (n - 2)<*- 3 ) 

< 

< C + 2d(n d - 2 + (n- 2) d ~ 2 + ... + e d ' 2 ) + 4d(n d " 3 + (n - 2) d ~ 3 + ... + e d ~ 2 ) 
<C + d{{n + l) d ~ 2 + n d ~ 2 + (n- l) d ' 2 + ... + l d ~ 2 ) + 4dn d ' 2 

rn+2 

<C + dx j t d ~ 2 dt + Adn d ~ 2 

= C+ -^—{n + 2)"" 1 - + 4dn d ~ 2 

d — 1 a — 1 

a — 1 

Therefore, 

r(n,d) < -^—n d - 1 + 0(n d - 2 ) d = 4,5,... (4) 
The following theorem summarizes our main results. 



12 



Theorem 3.1 For n > 1 and d>A, r(n,d) < ■£^n d ~ 1 + 0(n d_2 ). 
Specially for d = 4, f n 3 - f < r(n, 4) < |n 3 + 0(n 2 ) . 



4 Discussions 

In this paper we give an algorithm searching monotone ci-dimensional (d > 4) 
arrays A n> d, which requires at most -^n d ~ 1 + 0(n d ~ 2 ) comparisons. For d = 4, 
it is optimal up to the lower order terms. 

For d = 5, let D(n,5) = {(ii, i 2 , h, u, h) £ N 5 |*i + ^2 + h + U + ^5 = 
L|(n + l)J}U{(ii, i 2 ,i 3 ,i4,i5) e [nf\ii + i2 + iz + U + i$ = L|(n + 1)J + 1}, then 
\D(n, 5)| can be calculated to be ^n 4 + 0(n 3 ), which is the best lower bound 
on r(n, 5) currently known. However, applying the techniques in this paper, 
a i^-n 4 + 0(n 3 ) search algorithm for A n $ hasn't been found (our algorithm 
requires |n 4 + 0(n 3 ) comparisons in the worst case). So it may be interesting 
to tighten the bounds for d > 4. 
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